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. Abstract. For hamiltonian lattice gauge theory, we introduce the matrix product anzats inspired 

from density matrix renormalization group. In this method, wavefunction of the target state is 
assumed to be a product of finite matrices. As a result, the energy becomes a simple function 
of the matrices, which can be evaluated using a computer. The minimum of the energy function 
corresponds to the vacuum state. We show that the 5=1/2 Heisenberg chain model are well 
described with the ansatz. The method is also applied to the two-dimensional 5=1/2 Heisenberg 
and U(l) plaquette chain models. 
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f^. , RHIC experiments have started to test the fundamental properties of quantum chro- 

modynamics (QCD). The importance of first-principle analysis of QCD has increased 
largely in the context of color confinement. There is the expectation that quarks and glu- 
ons deconfine in extreme conditions such as heavy ion collision. However, lattice gauge 
theory at finite density has been stuck for a long time. The most ideal treatment of non- 
equilibrium quantum physics is to trace time-evolution of quantum states based on the 
Schrodinger equation. 

Density Matrix Renormalization Group (DMRG) is the variational method that gives 
the most accurate results in one-dimensional quantum systems LI] - In these days, DMRG 
has been used as a standard method to complement quantum Monte Carlo and suc- 
cessful in solving one-dimensional (zero and finite temperature) and two-dimensional 
(zero temperature) quantum systems. Application of DMRG to elemanrary particle and 
molecular physics has started some years ago. Recent interesting progress of DMRG is 
its application to non-equilibrium quantum physics and quantum information theory J2]. 

The first application of DMRG to particle physics was the massive Schwinger model 
with the term H3U . There is an old prediction by S. Coleman that quarks deconfine at 
6 = Tt [4]. The model has not been analyzed accurately with Monte Carlo because of the 
sign problem. On the other hand, DMRG has been successful in describing the details 
of the phase transition with large lattices because DMRG is free from the sign problem. 

The second application has been given by the author o]. It is a preliminary work 
for study of gauge theory. In bosonic lattice systems, each site has infinite degrees of 
freedom and therefore hamiltonian is infinite dimensional differently from spin and 
fermion systems. It is not evident whether DMRG truncation works for bosonic degrees 
of freedom. DMRG needs to be tested in a simpler bosonic model before going to gauge 
theory. In Ref. |5J], DMRG has been applied to a (l+l)-dimensional A0 4 model. The 
DMRG result for the critical exponent /3 = 0. 1264 ± 0.0073 is consistent with the exact 
one/3 = 1/8 = 0.125. 



TABLE 1. Numerical results for ground-state energy 
per site in the 5=1/2 Heisenberg chain model. The 
exact values have been obtained with the Bethe ansatz. 
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Matrix product variational method is a result of large simplification of DMRG |0]. 
In this method, the energy function has a simple form and easy to evaluate. However, 
the advantage is lost if a constraint is imposed to variational space directly. The Gauss 
law in gauge theory is one of the most important examples of constraints. To avoid this 
difficulty, the author has developed a method to introduce constraints with undetermined 
multipliers in Hamiltonian 0]. Also, the author has generalized the matrix product 
variational method so that it can be applied to higher-dimensional general systems Hi. 

We introduce matrix product states according to Ref. |6|]. With the knowledge from 
DMRG, wavefunction is represented as a product of matrices. 
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where periodicity is assumed. A [s] has the following normalization condition. 

M K 

L £ A * M «j8 A M a'jS = 8 aa > . (2) 
P=ls=l 

Each of the matrices A [s] can be parameterized with an appropriate number of indepen- 
dent variables |7]. It is expected that better results is obtained for larger M. Let us apply 
the ansatz to the 5=1/2 Heisenberg chain model with periodic boundary conditions 

H = t( S i S Ui+Spli+m + i)- (3) 

;=1 

Periodicity simplifies calculation of the energy function 

E[A] ~^W~\ tr(H) ' W 

where S a = Y, s .As\S a l \s')A*[s] ®A[s'] and 1 =E. s A*[s] ®A [s] . The minimum of the energy 
function corresponds to the ground state. In actual numerical calculation, the matrix 1 is 
diagonalized to simplify calculation of the powers of 1. 

Table[l]shows numerical results for ground-state energy per site, which are compared 
with the exact ones. M is the size of the matrices A [s\. When the lattice size L is small, 
convergence is poor even for large M. On the other hand, when L is large, the numerical 
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FIGURE 1. Reduction of two dimensional lattice to one dimensional 
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FIGURE 2. Energy per bond in the two-dimensional 5=1/2 Heisenberg model. 

results approaches the exact one as M becomes larger. In the best case, error is less than 
1%. The lattice size dependence of energy is small when L > 100. 

In the same way, the two-dimensional 5=1/2 Heisenberg model is analyzed. To 
use the matrix product ansatz, one-dimensional structure needs to be found on the two- 
dimensional lattice (see Fig. [T]). The tube and ring sizes are denoted as L\ and L2, 
respectively. Figure El plots energy per bond as a function of I/Z4 for a very large ring 
with L2 = 10000. The circle for Li = 1 is the exact result from the Bethe ansatz. The 
other three circles have been calculated using the matrix product ansatz. We are going 
to make a comparison in the thermodynamic limit L\ — > °° by extrapolating the obtained 
points to the limit. In Tabled energy per bond is compared among various results in the 
thermodynamic result. "Lattice size" means the largest lattice size used for calculation. 

The method is applied to a U(l) plaquette chain model. The hamiltonian is 

// = ££ / 2 _ x £ ([/p + [/ f )5 

/ p 

where x = l/g 4 H 1 311 . Figure|3]plots energy as a function of the parameter x and compares 
it with the results of Ref. [14]. In this calculation, the Gauss law V ■ E = has not been 
imposed on the variational space. Since vacuum wavefunction is available, we can check 
gauge invariance of the obtained vacuum state by calculating the vacuum expectation 
value of the electric field. The obtained result ~ 10~ 3 is small compared to energy, 



TABLE 2. Energy per bond of the two-dimensional 5=1/2 Heisenberg 
model in the thermodynamic limit. 
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FIGURE 3. Energy per site vs x in the U(l) hamiltonian lattice gauge model. 



which shows approximate gauge invariance of the vacuum state. Currently, calculation 
with the Gauss law constraint is being conducted. Further precise analysis with larger 
lattices will be given elsewhere to refine the results shown in this presentation. 

The numerical calculations were carried on the RIKEN RSCC system. This work has 
been partially supported by RIKEN BNL. 
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